Large thermomagnetic effects in weakly disordered Heisenberg chains 
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The interplay of different scattering mechanisms can lead to novel effects in transport. We show 
theoretically that the interplay of weak impurity and Umklapp scattering in spin- 1/2 chains leads to 
a pronounced dip in the magnetic field dependence of the thermal conductivity k at a magnetic field 
B ~ T. In sufficiently clean samples, the reduction of the magnetic contribution to heat transport 
can easily become larger than 50% and the effect is predicted to exist even in samples with a large 
exchange coupling, J > B, where the field-induced magnetization is small. Qualitatively, our theory 
might explain dips at B ~ T observed in recent heat transport measurements on copper pyrazine 
dinitrate, but a fully quantitative description is not possible within our model. 
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Some of the most fascinating manifestations of quan- 
tum many-body physics in one-dimension (ID) can be 
found in spin-i chain systems^. In particular, the 
spin- ^ Heisenberg model with antiferromagnetic nearest 
neighbor exchange interactions is one of the most exten- 
sively studied paradigms. It provides a remarkably non- 
trivial example of an exactly solvable mode£, allowing a 
detailed analysis of its thermodynamic properties. While 
integrability is special to this particular model, its low en- 
ergy properties are generic: they do not differ essentially 
from low energy properties of other (non integrable) spin 
chain models, with more general finite range interactions. 
In essence, these systems support Fermionic elementary 
excitations - spinons - which carry spin and no charge. 
Their kinetic energy and interactions are dictated by the 
exchange couplings in the chain, and their Fermi mo- 
mentum can in principle be tuned by a magnetic field B, 
which plays the role of a chemical potential. 

Thermodynamic properties of low-dimensional spin 
systems and especially of the Heisenberg model are gen- 
erally very well understood allowing for a quantitative 
description of a broad range of experiments. In compar- 
ison, the heat and spin transport 4-5 in spin-systems is 
considerably more complicated. For example, the heat 
conductivity k of the one-dimensional Heisenberg model 
is infinite at finite temperatures as the heat current op- 
erator is a conserved quantity for this idealized model. 
In real materials, the effects of disorder, phonon cou- 
pling and spin interactions not captured by the integrable 
Heisenberg model render the conductivity finite. Never- 
theless, it remains often very large& as long as disorder is 
weak. 

In this paper we study how the strong interactions of 
the Heisenberg model affect the heat transport in the 
system in the presence of weak disorder beyond the well 
known effect that disorder is strongly renormalized by 
the interactions 7 -. Our study is directly motivated by re- 
cent experiments 5 - in the spin- 1/2 chain compound cop- 
per pyrazine dinitrate Cu(C4H 4 N 2 )(N03)2 (CuPzN) re- 
produced in Fig. [TJ In this system a small exchange 



coupling, J/ks ~ 10.3 K, allows to polarize the sys- 
tem with moderate magnetic fields B. While both the 
spins and the phonons contribute to the heat transport 
at low temperatures T with similar strength, the field 
dependence can be used to separate the two effects. The 
total heat conductivity k can be split into a contribu- 
tion arising purely from phonons and a magnetic part, 
k(B,T) = K p h(T) + K mag (T, B). Then, one can use the 
fact that K p h(T) is practically independent of B to ex- 
tract 

An mag {T,B) = k(T, B) - k(T, B = 0) 

= K mae (T,B)-K maK (T,B = 0). (1) 

In CuPzN AK mag (T, B) shows a pronounced dip as a 
function of magnetic field for small B and T. Inter- 
estingly, the position of this dip scales linearly with the 
temperature, g^BB min w SksT, pointing to a simple un- 
derlying mechanism. A dip in n ma ,g(B) has been observed 
in numerical simulations 8 of classical spin-chains coupled 
to phonons but this dip occurs at B ~ J and does not 
scale with T. 

Several processes can contribute to K mag , the field de- 
pendent part of k. First, there is a positive contribution 
from heat transported by the spin-chains. Second, the 
heat conduction of the phonons is reduced when phonons 
scatter off spin fluctuations. Third, there is a contribu- 
tion from spin-phonon drag which is usually positive and 
can also become very large 9 . The positive sign of K mag 
in CuPzN suggests that the first and possibly the third 
mechanism are dominating. 

Within a fermionic interpretation of the spin excita- 
tions, the magnetic field enters as a chemical potential 
while T determines the broadening of the Fermi sur- 
face. Therefore a likely interpretation of the experiment 
is that the characteristic dip arises when the Fermi sur- 
face moves away from the momentum kp = n /2a, cor- 
responding to a half-filled system with lattice spacing a. 
The commensurate filling for B = is very special as 
it allows for Umklapp scattering: two spinon excitations 
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FIG. 1: (Color online). Upper figure: Field dependent 
thermal conductivity k(B,T) — k(0, T) in the spin-1/2 com- 
pound CuPzN taken from Ref. [5|]. The downturn at large 
fields arises as the spinon band is considerably depleted. For 
T — the magnetization saturates at B « 15 T. Lower figure: 
The normalized magnetic part of the thermal conductivity, 
Kma g (-B,T)/K mag (0,T), as a function of h = fi B gB/(k B T) 
shows a dip at h ~ 3. Note that there might be a consider- 
able error in the size of the dip as K ma g(S = 0, T) could not 
be measured directly due to a large phonon background, see 
text and Ref. [f| for details. 

can be transferred from the left to the right Fermi surface 
and the excess momentum Akp — 2n/a can be absorbed 
by the underlying lattice. As this scattering mechanism 
is only effective in the vicinity of the commensurate point, 
it is exponentially suppressed for gfisB > knT . We will 
show that besides the Umklapp scattering it is necessary 
to include the effects of impurities to get structures at 
B ~ T; in the absence of disorder, the presence of certain 
approximate conservation laws 10 prohibits a relaxation of 
the heat current by the leading Umklapp process [see §5§ 
below] alone. 

Model and methods: We investigate the one- 
dimensional Heisenberg model in the presence of a mag- 
netic field B and weakly disordered exchange couplings, 
8J t < J 

H = -^(J + SJjSi-Si+t-gUBB^Sz. (2) 

i i 

In the following, we will be interested in the limit of small 
magnetic fields and temperatures, B, T <C J. More pre- 
cisely, all calculations are done only to leading order in 
I/ln[(T, B)/J]. The interactions lead to a strong renor- 



malization of disorder—, but we assume that the temper- 
atures are sufficiently high, such that the renormalized 
disorder remains weak, 5 J <C V JT. 

For B,T,5Ji <C J one can use the powerful techniques 
of bosonization to describe the low-energy properties of 
the system. It is useful to split the effective low-energy 
Hamiltonian into three pieces 

H = Hll + Hu + H c ii S (3) 

Hll = v J^-(K{d x 9f + ^{d x ^ (4) 

H v = j^-^ J dx[e lAkx e l4 * + h.c] (5) 

H D = -!- / dxr 1 (x)[ie* 2 ' t> + h.c] (6) 

27TO J 

where (using the notation of 011(3]) d x <t> denotes fluctua- 
tions of the magnetization in the z direction, 9 is the con- 
jugate variable with [4>{x), d x 9{x')] = iir5{x — x'), and we 
use units where ks = 1, fi = f • For the spin-rotationally 
invariant Heisenberg model, the velocity and the Lut- 
tinger parameters at the low-energy fixed point are given 
by v — \ Ja and K = 1/2, respectively. The Umklapp 
term Hjj describes the scattering of spinons from one 
Fermi point to the other. At a finite magnetization, the 
Fermi momentum kp = ^■(l + 2(5' z )) deviates from 7r/2a 
and therefore the excess momentum Afc = 4fcp — — = 

1 a 

An(S z )/a cannot be absorbed by the crystalline lattice. 
For Ak — the Umklapp scattering is a marginally irrel- 
evant operator whose strength decreases logarithmically 
with T (see below) while for vAk 3> ksT it is exponen- 
tially suppressed in a clean system. The effects of disor- 
der, or more precisely from components of the disorder 
potential oscillating with momentum 2kp, is described by 
Ho which models the scattering from one Fermi point to 
the other with rj{x) ~ 8J{x). Here we assume uncorre- 
cted disorder with (r](x)ri(x')) = D,n s 8{x — x 1 ). 

To calculate the heat conductivit y w e use the so-called 
memory matrix formalismii as in [10(. Within this ap- 
proach one calculates a matrix of relaxation rates for a 
given set of modes. As has been discussed in detail in 
[13], in general this method allows to calculate a lower 
bound to the conductivity. The formalism gives precise 
results as long as the relevant slow modes are included in 
the calculation. For the present system, the essential step 
is to realiz o 10 i 13 that in the absence of disorder the oper- 
ator Q = Jh + v fjfJs is conserved, [Q, H LL + Hu] = 0. 
Here, Jh = v 2 J dxd x 8d x <fi is the heat current associ- 
ated with Hll and J s — vK J 89/ ir is the spin-current. 
This can be seen by realizing that up to the prefactor 
v 2 the heat current can be identified with the momen- 
tum operator, the generator of translations. The Umk- 
lapp term describes a process where a momentum AA: is 
generated and therefore its commutator with Jh is pro- 
portional to v 2 Ak. Similarly, the spin current is changed 
by — Av as two spinons with velocity v are scattered into 
states with velocity —v. Therefore the linear combina- 
tion Q = Jh + VjjtJ s remains unaffected by Umklapp 



processes. In terms of the original variables, Q can be 
identified with J 2 ^\ S.j(S i+1 x S i+2 )i the heat current 
operator for B = which commutes with the integrable 
Heisenberg model ([2} in the absence of disorder, SJi = 
(when longer range interactions or interchain coupling 
break integrability, the lifetime of Q nevertheless remains 
exponentially large in a clean system, see 0IH). We 
therefore set up the memory matrix formalism in oper- 
ator space spanned up by the two relevant currents J B 
and J s . 

The decay rates of the currents are determine d 1 ^ 11 
from a 2 x 2 matrix M of correlation functions (the 'mem- 
ory matrix'), 
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(7) 



where (A; B) u denotes a retarded correlation function of 
A and B evaluated at the real frequency u and J\ — 
Jh, J2 = Js are the two relevant operators. Within the 
assumptions of our paper, we can treat both Umklapp 
and disorder perturbatively and therefore it is sufficient 
to evaluate all expectation values with respect to Hll 
with K = 1/2 as the derivative dtJ% are already linear in 
the perturbations Hjj and H^is- 

The heat conductivity per spin chain is obtained from 
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T M SS M HH 

is the susceptibility of 



where \h = {Jh;Jh)u=o 
the heat current. 

Separating the contributions from Umklapp and dis- 
order, M = Mjj + Mdis, we find using straightforward 
perturbation theory 



M v = 

?u = 
Mdis = 




g(r)MAfc) : 

8tt 2 

vnDdis ( T 
8a i 4 



n' B (vAk/2) 



(9) 
(10) 
(11) 



where n' B (uS) is the derivative of the Bose function, 
n B (u) = l/(e w / T - 1). The T dependence of g(T) ~ 
irv/\n(J/T) (see below) takes into account that the Umk- 
lapp scattering is marginally irrelevant with respect to 
the clean fixed point. 

Results: The rather simple equations (|51TTTj) describe 
the complex interplay of disorder and Umklapp scatter- 
ing. First, in the absence of disorder, the heat conductiv- 
ity is infinite 14 as Mjj has an eigenvalue reflecting the 
conservation of Q described above. Second, for vanish- 
ing magnetic field, Afc — 0, Umklapp scattering plays no 
role and one obtains the well-known result K mag /T ~ T. 
The mean free path decreases linearly in T as the disor- 
der is strongly renormalized by the interactions 7 . Third, 




FIG. 2: (Color online). Field dependence of the of the nor- 
malized thermal conductivity as a function of the rescaled 
field, h = g/iBB/ksT. For moderate disorder a(T), Eq. (fT2J, 
10 > ot(T) > 1, a pronounced dip in K mag (-B) is predicted for 
gfJ,BB ~ 4ifiT. For even smaller values of a{T) the dip gets 
broader and deeper, approximating the asymptotic behavior 
(|14[) (solid black line), up to larger and larger fields. 



for finite magnetic field and T — * 0, Umklapp scatter- 
ing is exponentially suppressed, IV ~ e -vAk/2T^ ag ^ e 
Fermi energy has shifted away from commensurate fill- 
ing and one finds K m&g (B ^> T) = K ma . s {B = 0). Here 
we have neglected the - formally subleading - effect that 
the Luttinger liquid parameter K and therefore also the 
renormalization of the disorder potential depend on B. 

Upon increasing the magnetic field B, the Fermi sur- 
face is shifted and Afc increases approximately linearly 
in B, Afc AirxB/a with 1 ^ x ~ 9^b/{t^ 2 J)- As argued 
above, for Afc = the Umklapp scattering does not in- 
fluence transport. Therefore, by raising Afc the effect 
of Umklapp scattering is switched on proportionally to 
(Afc) 2 . But upon increasing Afc further, Umklapp scat- 
tering is switched off for uAfc > T or B > T. The net 
result is a pronounced dip in « mag (i?), see Fig. [2] 

In the scaling limit of weak disorder and 
1/ ln[(B, T)/ J] <C 1 the normalized conductivity 



l magV 



(B,T)/K ma , s (B = 0,T) is only a function of the 
h = [isgB / ksT and a dimensionlcss 



scaling variable 
variable 



a(T) = 



D dis v 2 



(k B T) 2 g(T) 2 a 9irg(T) 2 K mag (B = 0, T) 



(12) 



which parameterizes the relative strength of (renormal- 
ized) disorder and Umklapp scattering (g ~ 1/ ln[ J/T] is 
defined below). In these variables, we obtain 



TT 3 a(T)-2TT 2 h 2 n' B (h) 



«ma g (0, T) n 3 a{T) - (2ir 2 + Ah 2 )h 2 n' B {h) 



(13) 



with n B (h) = l/(e h - 1). As shown in Fig. [2j the field 
dependence of the thermal conductivity is predicted to 
show a pronounced dip at B ~ T. For small a(T) <C 
h 2 n' B (h), i.e. weak disorder and not too strong fields, 
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one finds 

KmagjB, T) ^ 1 , . 

«mag(0,T) ~ 1 + ^ ' 

see Fig. [5J This implies a strong reduction of K mag of 
order 1 for [IbB ~ fc^T as long as the renormalized dis- 
order is sufficiently weak, a(T) <C 1. This is the main 
result of the paper. For large fields or stronger effective 
disorder one obtains a small suppression of K mag , 

K ma8 (B,T) ^ 1 -4r^(fe)fe 4 
(0, T) ir 3 a(T) 

giving rise to a minimum at h ~ 3.8 of size 0.63/a(T). 

To allow for a quantitative comparison to the experi- 
ment of Ref. [f|, one needs to estimate a(T). For this, 
both the strength of impurity scattering and the strength 
of the renormalized Umklapp scattering g(T) have to be 
determined. Fortunately, for simple Heisenberg chains 
the latter is known analytically from the Bethe Ansatz. 
Translating our notations to those used in [l5j], we ob- 
tain g(T) = g(T)Tr 2 Ja/2 and g is obtained by solving the 
equation^ 

i + ^f>=I»p^) (16, 

where 7 = 0.5772... is the Euler constant. Within the 
precision of our calculation, g ss l/ln(J/T), but we 
use the more precise formula (|16p which includes sub- 
leading corrections below. The disorder strength for 
the sample of CuPzN measured in 0] can in princi- 
ple be obtained from n mag (B = 0,T). Unfortunately, 
a large phonon background prohibits a direct measure- 
ment of this quantity but a crude estimate, K mag w 
3.5 T 2 Wm _1 K~ 3 , can be obtained from the behavior of 
K at large fields (see @ for details). Using this esti- 
mate, we find for the heat conductivity per spin chain 
K mag w 2.1 l(r 18 T 2 WmK" 3 . For the four lowest tem- 
peratures, T = 0.37,0.66,0.96, 1.48K, shown in FigHJ 
we obtain from Eq. (JT2J) a(T) « 0.52,0.12,0.049,0.016, 
respectively. 

These estimates allow a quantitative comparision of 
theory and experiment. There are two main discrepan- 
cies between theory and experiment which can be seen 
from a direct comparison of Fig. [T] and Fig. [2j First, 
there is a discrepancy in the position of the minimum 
(located at h w 3 in the experiment and at h « 4 within 
the theory) and second, the size of the dip of the order 
of 10% is much smaller than the predicted reduction of 
more than 50% - or the estimate for a appears to be 
almost two orders of magnitude too small. What can be 
the origin of the clear discrepancy? First, one should note 
that for the temperatures and magnetic fields shown in 
Fig. ED both subleading effects of order \n{B/T) / ln( J/T) 
or ln(ln( J/T))/ ln( J/T) and band-curvature effects (the 
overall downturn of K mag in large fields) neglected in our 



calculation can become important. For example, g 2 cal- 
culated to leading order is for J/T = 30 a factor 2.5 
larger than the value obtained from Eq. Q16p . 

More importantly, we believe that our model © does 
not capture all aspects of the physics in the CuPzN 
samples correctly. Especially, modeling the disorder by 
Eq. ([5]) might not be appropriate. This was also the 
conclusion of Ref. Q from an analysis of the heat con- 
ductivity at large fields B ~ 15 T in the quantum critical 
regime where the magnetization is close to saturation. 
Indeed, for other types of disorder the matrix pip will 
have a different structure which will affect the quanti- 
tative predictions while the qualitative picture will re- 
main unmodified. For example, it might be necessary to 
take the interplay of forward scattering from impurities 
and interactions into account. Forward scattering affects 
transport at B — only very weakly but can reduce the 
Umklapp dip in ft mag considerably as the suppression of 
«mag at larger fields relies on momentum conservation. A 
more realistic modeling of disorder should also account 
for the possibility that defects cut the one-dimensional 
chains in long pieces^. In such a situation, one has also 
to model how phonons (or weak inter-chain interactions) 
couple heat into and out off such long chain segment a 5 ' 16 . 

Conclusions: Our theoretical calculations show that 
the heat conductivity of weakly disordered spin chains 
is very sensitive to moderate magnetic fields B ~ T. A 
pronounced dip in the field dependence of k arises from 
the shift of the Fermi surface of the spinons induced by 
the magnetic field. The effect of Umklapp scattering on 
the heat conductivity turns out to be strongest when the 
Fermi surface is shifted from the commensurate position 
at B = by an amount of the order of its thermal broad- 
ening. 

It is interesting to note that, according to our theory, 
this effect should be observable for a wide range of pa- 
rameters including spin-chains which have in contrast 
to CuPzN - large exchange couplings of several hundred 
Kelvin. Due to the strong B dependence, it should be 
possible to identify dips in K mag even in the presence 
of a large phonon background. However, for such sys- 
tems, the effective disorder has to be sufficiently small, 
a(T) < 10. To obtain an effective disorder strength of 
the order of 1 at B ~ T, typical fluctuations of the ex- 
change coupling have to be of the order of B (for this 
crude estimate we used -Ddis ~ {SJ) 2 a and neglected 
logarithmic renormalizations) 6 J / J < (J-bB/ J. Further- 
more, Afc ~ hbB/(Jo) is also small for large J and a 
necessary condition for the quantitative validity of our 
calculations is that there is no substantial forward scat- 
tering on the associate length scale 1/Afc ~ aJ / '(ij,bB). 
In systems with large J and strong phonon scattering, 
one has also to take into account^ that the sound veloc- 
ity c is smaller than the spinon- velocity v. Therefore it 
may happen that the position of the dip in K mag moves to 
lower values, h ~ c/v, as the relevant energy scaled for 
phonon-assisted Umklapp scattering is cAfc rather than 
vAk. 
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